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Abstract 

The extrapolation of small-cluster exact-diagonalization calculations is used 
to study the influence of doping on valence transitions in the spinless Falicov- 
Kimball model at nonzero temperatures. Two types of doping are examined, 
and namely, the substitution of rare-earth ions by non-magnetic ions that in- 
troduce (i) one or (ii) none additional electron (per non-magnetic ion) into the 
conduction band. It is found that the first type of substitution increases the 
average /-state occupancy of rare-earth ions, whereas the second type of substi- 
tution has the opposite effect. The results obtained are used to describe valence 
transition behavior of samarium in the hexaboride solid solutions Sm\- X M X B§ 
(M = y 3+ , La 3+ , SV 2+ , Yb 2+ ) and a very good agreement of theoretical and 
experimental results is found. 
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In the past decade, a considerable amount of effort has been devoted to un- 
derstanding of the multitude of anomalous physical properties of rare-earth and 
transition-metal compounds. Recent theoretical works based on exact numerical 
and analytical calculations showed that many of these anomalous features, e.g., 
mixed valence phenomena, metal- insulator transitions, etc., can be described very 
well within different versions of the Falicov-Kimball model (FKM) P3 |2l E]. In 
particular, it has been found that the spinless FKM, can describe the discontinuous 
valence and insulator- metal transitions jUE] (induced by external pressure) observed 
experimentally in some rare-earth sulfides, halides and borides jB] . A nice correspon- 
dence between theoretical and experimental results has been found, for example, for 
SitiBq |Hl E] , where the alternative explanation of formation and disappearance 
of the energy gap in the electronic spectrum of SttlBq has been presented. For StuBq 
a very good agreement of experimental and theoretical results has been obtained also 
for nonzero temperatures for both the temperature dependence of the specific heat 
as well as magnetic susceptibility jTUj. In addition, it was found that the spinless 
FKM can describe the anomalous temperature increase of electrical conductivity a 
in StuBq (due to the valence transition) ^T] as well as the low temperature behavior 
of a (due to the appearance of the in-gap structure) in this material E]. Re- 
cently, the FKM has been successfully used to describe the inelastic light scattering 
in the StuBq systems [T3|. These results indicate that the spinless FKM is a quite 
good microscopic model that can describe, at least qualitatively, many of anomalous 
features of rare-earth compounds. From this point of view it seems to be interesting 
to ask if this simple model can describe also another kind of experiments performed 
often in these materials, and namely, the influence of doping on the valence transi- 
tion behavior of systems. For SmB 6 it has been found, for example fUE3> that the 
substitution of Sm by non-magnetic divalent ions (e.g., Sr 2+ ,Yb 2+ ) increases the 
average samarium valence (the average occupancy of / orbitals decreases), whereas 
the substitution of Sm by non-magnetic trivalent ions (e.g., Y 3+ ,La 3+ ) produces 
the opposite effect. In this paper we try to describe the influence of both types of 
substitutions on the average samarium valence within the spinless FKM. 

The spinless FKM is based on the coexistence of two different types of electronic 
states in given materials: localized, highly correlated ioniclike states and extended, 
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uncorrelated, Bloch-like states. It is generally accepted that valence and insulator- 
metal transitions result from a change in the occupation numbers of these electronic 
states, which remain themselves basically unchanged in their character. The Hamil- 
tonian of the model can be written as the sum of three terms: 

# = E MM + u £ ftfidfd, + E f '£f+f i , (1) 

ij i i 

where /+, /j are the creation and annihilation operators for an electron in the lo- 
calized / state of a rare-earth atom at lattice site % with binding energy Ef and df, 
di are the creation and annihilation operators of the itinerant spinless electrons in 
the ci-band Wannier state at site i. The first term of (1) is the kinetic energy corre- 
sponding to quantum-mechanical hopping of the itinerant d electrons between sites 
i and j. These intersite hopping transitions are described by the matrix elements 
tij, which are —t if % and j are the nearest neighbors and zero otherwise (in the 
following all parameters are measured in units of t). The second term represents the 
on-site Coulomb interaction between the <i-band electrons (with the total number 
Nd = Y^idfdi) and the localized / electrons (with the total number Nf = Y^iftfi)- 
The third term stands for the localized / electrons whose sharp energy level is Ef. 
The Hamiltonian of the FKM doped by divalent (e.g., Sr 2+ ) or trivalent (e.g., Y 3+ ) 
ions has the same form. From the numerical point of view there is only one differ- 
ence, and namely, the summation over all lattice sites (in the second and third term) 
should be replaced by summation over the samarium positions. 

Since in the spinless version of the FKM without hybridization the /-electron 
occupation number f^fi of each site i commutes with the Hamiltonian (1), the f- 
electron occupation number is a good quantum number, taking only two values: 
Wi = 1 or 0, according to whether or not the site i is occupied by the localized / 
electron. Then the Hamiltonian (1) can be written as 

H = Y, h v d tdj + E f Y,w u (2) 

ij i 

where hij{w) = + UwiSij. 

Thus for a given /-electron configuration w = {wi,w 2 ■ ■ ■ wl} defined on the 
one or two-dimensional lattice of L sites with periodic boundary conditions, the 
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Hamiltonian (2) is the second-quantized version of the single-particle Hamiltonian 
h(w) = T + UW, where T is an L-square matrix with elements tij and W is the 
diagonal matrix with elements Wi. The behavior of this model at nonzero temper- 
atures is studied using small-cluster exact-diagonalization calculations. In order to 
compensate partially for the small size of clusters, thermal properties of the system 
are investigated via the grand canonical ensemble. 

The grand canonical partition function of the FKM ^H] can be expressed directly 
in terms of the eigenvalues of the single-particle operator h, which depend on the 
/-electron configuration specified by w = {wi}, i.e., 

S = £ e-^t-riW f[[l+ e-^ )/T \. (3) 

w i=l 

The temperature dependences (r = ksT) of the thermodynamic quantities are 
calculated from this expression by using the standard statistical mechanics procedure. 
For example, the total number of /-electrons Nf and the internal energy (S) as 
functions of chemical potential \i are given by 



where /3 —1/r. 

In the next, we consider (i) the substitution of rare-earth ions (e.g., Sm) by 
non-magnetic trivalent ions (e.g., Y 3+ ) which introduce conduction electrons into 
the (^-conduction band (one electron per dopant) and (ii) the substitution of rare- 
earth ions by non-magnetic divalent ions (e.g., Sr 2+ ) which play a dilution role and 
reduce the number of conduction electrons in the ^-conduction band (no additional 
electrons are introduced to the system). In the first case the total number of electrons 
N = Nf + Nd is equal to the total number of lattice sites L (the case of one electron 
per rare-earth atom is considered) while in the second case N = L — X, where 
X is the total number of dopants. Thus the first step of numerical calculations 
is to determine the chemical potential corresponding to N = L and N = L — X , 
respectively. Afterwards rif, £ and other thermodynamic quantities can be calculated 
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numerically. 

It should be noted that the summation over all / electron distributions in (3) 
should be replaced for the doped FKM by summation over all possible /-electron and 
dopant distributions, thereby strongly limiting numerical computations. Although 
the number of configurations can be reduced considerably by the use of symmetries 
of H, there is still a limit (L ~ 20) on the size of clusters that can be studied 
with this method. However, we will show later that due to small sensitivity of the 
thermodynamic quantities of the FKM on L, already such small clusters can describe 
the thermodynamics of the model very well. 

To describe effects of doping by divalent and trivalent ions on the average /-state 
occupancy of rare-earth ions we have performed an exhaustive study of the model for 
a wide range of model parameters U, r, Ef. The typical dependences of the average 
/-orbital occupancy nf = Nf/(L — X) on the dopant concentration x = X/L are 
displayed in Fig. 1 for Ef = 0, U = 0.5 and several different temperatures r. 




Figure 1: The average occupancy of / orbitals as a function of x for doping by diva- 
lent and trivalent ions calculated for L = 12, 16, 20 and several different temperatures 
r. Curve a, t = 0.1; curve b, r = 0.3; curve c, r = 0.6; curve d, r = 1. Inset: The 
average samarium valence in the Smi- X M X BQ systems (M = Y 3+ , La 3+ , Sr 2+ , Yb 2+ ) 
measured at 300 K (Ref. 14). 
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To determine the finite-size effects all concentration dependences have been cal- 
culated for several cluster sizes (L = 12, 16, 20) and displayed on the common linear 
scale. We have found that for all examined temperatures r the finite-size effects are 
negligible (over the whole range of x plotted) and thus these results can be used 
satisfactorily to represent the behavior of macroscopic systems. It is seen that the 
substitution of rare-earth ions by trivalent and divalent ions produces fully different 
effects. While in the first case the average occupancy of /-orbitals increases with 
increasing dopant concentration x, in the second case nj decreases. In all examined 
cases the valence transition realizes continuously. Comparing these results with ex- 
perimental measurements of the average samarium valence v (see inset in Fig. 1) in 
the Sm 1 _ x M x B 6 systems HHE3 ( M = y3+ > La3+ > Sr2+ , Yb 2+ ) one can find a good 
qualitative agreement between the theoretical and experimental results. Indeed, the 
substitution of Sm by non-magnetic divalent ions (Sr 2+ ,Yb 2+ ) increases the aver- 
age samarium valence (the average occupancy of / orbitals rif = 3 — v decreases), 
whereas the substitution of Sm by non-magnetic trivalent ions (Y 3+ , La 3+ ) produces 
the opposite effect. These results again indicate that the FKM, in spite of its rela- 
tive simplicity could, in principle, yield the correct physics for describing rare-earth 
compounds. Even, also small quantitative deviations between the experimental and 
theoretical results can be straightforwardly improved within the FKM by optimizing 
the rif(x) behavior with respect to r and Ef. This is illustrated in Fig. 2, where the in- 
dependence of the /-orbital occupancy is plotted for r = 1 (this relatively large value 
is chosen to ensure that r is much greater than the spacing of many-body energy- 
levels for both D = 1 and D = 2) and several different values of Ef. It is seen that 
for nonzero values of Ef a nice correspondence between the theoretical and experi- 
mental results can be achieved. Of course, one can take objection that real materials 
are three dimensional, while our results have been obtained for the one dimensional 
case. Unfortunately, due to limitations (mentioned above) on the size of clusters that 
can be treated exactly by small-cluster exact-diagonalization calculations (L = 20) 
we are not able to obtain any exact results in three dimensions. However, in two 
dimensions there are several finite clusters with L = 10, 16, 18, 20 that are accessi- 
ble to exact-diagonalization calculations ^B] and can be used to test the stability 
of obtained results with increasing dimension. The numerical calculations that we 
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Figure 2: Comparison of the one-dimensional theoretical (dashed lines) and ex- 
perimental results for the average Sm valence v = 3 — n/ as a function of x. The 
theoretical results have been calculated for U = 0.5, r = 1 and several different values 
of the /-level position Ef. The finite-size effects are negligible, on the linear scale it 
is not possible on the drawing to distinguish behaviors obtained for L = 12, 16 and 
20 over the whole range of x plotted. The experimental results have been obtained 
for T=300 K. 

have performed in two dimensions (for the same set of model parameters as in D=l 
and the above mentioned cluster sizes) showed (see Fig. 3) that the /-occupancy 
of rare-earth orbitals depends only very weakly on the dimension of the system. 
Thus we can conclude that the spinless FKM, when treated exactly, can provide a 
reasonable theoretical description of valence transitions induced by doping in some 
rare-earth compounds, e.g., the samarium hexaboride solid solutions Sm\_ x M x B§ 
(M = Y 3+ , La 3+ , Sr 2+ , Yb 2+ ). Of course, this statement can not be extended auto- 
matically on all thermodynamics characteristics of the doped FKM. For example, to 
do quantitative comparison between the theoretical and experimental results for the 
specific heat and the magnetic susceptibility, the model should be first generalized 
by including the spin and orbital dynamics that description was too simplified in the 
Hamiltonian (1). The work in this direction is currently in progress. 

In summary, the extrapolation of small-cluster exact-diagonalization calculations 
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Figure 3: Comparison of the two-dimensional theoretical (dashed lines) and ex- 
perimental results for the average Sm valence v = 3 — n/ as a function of x. The 
theoretical results have been obtained on clusters with L = 10, 16, 20 for the same 
set of U, r and Ef values as in D = 1. Again, the finite-size effects are negligible. 

was used to study the influence of doping on valence transitions in the spinless FKM 
at nonzero temperatures. Two types of doping were examined, and namely, the 
substitution of rare-earth ions by non-magnetic ions that introduce (i) one or (ii) 
none additional electron (per non-magnetic ion) into the conduction band. It was 
found that the first type of substitution increases the average /-state occupancy of 
rare-earth ions, whereas the second type of substitution has the opposite effect. In 
all examined cases valence changes are continuous. The results obtained were used 
to describe valence transition behavior of samarium in the hexaboride solid solutions 
Smx- x M x Be (M = Y 3+ , La 3+ , Sr 2+ , Yh 2+ ) and a very good agreement of theoretical 
and experimental results was found. 
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